from Excitation_Carbon import *

cs = CS_Excitation(species_name = "C+5")


id = 0

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =367.0,
    A   = 2.446e-2,
    B   = -1.567e-2,
    C   = 1.010e-2,
    D   = 0.0,
    E   = 0.0,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+5_1s-2s")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =367.0,
    A   = 1.0406e-2,
    B   = 1.720e-2,
    C   = 4.370e-2,
    D   = 0.0,
    E   = 1.259e-1,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+5_1s-2p")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 435.0,
    A   = 4.866e-3,
    B   = -2.522e-3,
    C   = 3.158e-3,
    D   = -1.765e-3,
    E   = 0.0,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+5_1s-3s")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =435.0,
    A   = 4.248e-3,
    B   = 8.978e-3,
    C   = 4.082e-3,
    D   = 0.0,
    E   = 2.130e-2,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+5_1s-3p")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =435.0,
    A   = 3.874e-3,
    B   = -6.100e-3,
    C   = 3.909e-3,
    D   = 2.143e-3,
    E   = 0.0,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+5_1s-3d")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =459.0,
    A   = 1.814e-3,
    B   = -1.035e-3,
    C   = 1.515e-3,
    D   = -9.157e-4,
    E   = 0.0,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+5_1s-4s")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =459.0,
    A   = 1.106e-3,
    B   = 2.965e-3,
    C   = 3.247e-3,
    D   = 0.0,
    E   = 7.915e-3,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+5_1s-4p")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =470,
    A   = 6.869e-4,
    B   = 1.166e-3,
    C   = 2.120e-3,
    D   = 0.0,
    E   = 3.711e-3,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+5_1s-5p")



cs.plot_cs()
cs.save_cs()
cs.cal_total_cs()